Ionization of the Xenon Fluorides 
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The noble-gas atom xenon can bind fluorine atoms and form several stable compounds. We study 
the electronic structure of the xenon fluorides (XeF n , n = 2,4,6) by calculating their ionization 
spectra using a Green's function method, which allows one to treat many-body effects at a high 
level. Our focus is on the valence region and on the Xe4d core hole. We observe a sensitive 
dependence of the spectra on the number of fluorine ligands. Systematic line shifts are uncovered 
and explained. In the Xe5s and F2s inner- valence regimes, from XeF2 to XeF6, the usefulness of 
the one-particle picture of ionization is found to become progressively worse. Moreover, adding the 
electronegative fluorine ligands seems to enhance — not suppress — the Auger decay of a Xe4ci core 
hole. 



I. INTRODUCTION 

Noble-gas atoms are, in general, chemically rather in- 
ert. Nevertheless, since the heavy noble-gas species kryp- 
ton, xenon, and radon possess a comparatively small ion- 
ization potential, they can form molecules, at least to- 
gether with electronegative elements. The most stable 
and most widely investigated among these compounds 
are the xenon fluorides 1 XeF2, XeF4, and XeFg. 

The electronic structure of the xenon fluorides has been 
the subject of a number of both experimental and theo- 
retical investigations. In the 1970's, photoabsorption ex- 
periments, at photon energies between 50 and 160 eV 2,3 
and between 6 and 35 eV 4 , were carried out at DESY in 
Hamburg. Carroll et al. 5 employed photoelectron spec- 
troscopy to measure the chemical shifts of the F Is and 
Xe 3d levels in the xenon fluorides. More recently, using 
synchrotron radiation, Cutler et al. 6 obtained gas-phase 
photoelectron spectra with such a high resolution that 
fluorine-ligand field splittings on the Xe4d levels could 
be extracted. 

Photoelectron spectroscopy, in combination with ab 
initio studies, is an outstanding tool for characterizing 
the electronic structure of a molecule. However, at least 
in the case of the xenon fluorides, there exist — to our 
knowledge only few theoretical contributions to moti- 
vate and back-up experimental work. The pioneering pa- 
pers 7-9 , which are more than 20 years old, utilized self- 
consistent-field methods to calculate ionization poten- 
tials. The restriction to an effective one-particle model 
made it impossible to uncover complexities due to elec- 
tron correlation. Newer theoretical work 10-12 on xenon 
fluorides focused on electronic ground-state geometries 
and dissociation energies. It is therefore timely to im- 
prove on the early studies and to explore some of the 
many-body physics of electrons in XeF2, XeF4, and XeFg. 
This is the purpose of this paper. 

We proceed as follows. Section II describes the Green's 
function method we use to calculate ionization spectra. 
We discuss the molecular geometries we employed in our 
calculations, and we analyze the impact of relativistic 
effects on the ionization spectra of the xenon fluorides. 



The ionization spectra we calculated cover the valence 
regime and the Xe4d core line; they are presented in 
Sec. III. A detailed comparison between spectra that 
take electron correlation into account and spectra based 
on the Hartree-Fock model demonstrates the occurrence 
of many-body phenomena in almost the entire spectral 
regime we consider. These phenomena display an inter- 
esting dependence on the number of fluorine ligands. A 
summary and conclusions are given in Sec. IV. 



II. COMPUTATIONAL DETAILS 
A. Algebraic Diagrammatic Construction 

Green's functions are a fundamental tool of many-body 
theory. 13-16 They are well suited to calculating various 
properties of molecules. 17,18 For instance, the pole po- 
sitions of the one-particle propagator yield the vertical 
ionization potentials of a molecule. The residuum or 
pole strength of a pole of the one-particle propagator 
is a measure of the one-particle character of a specific 
ionized state, i.e., it is a measure how well the physical 
state is described within the one-hole configuration space 
deriving from the Hartree-Fock ground state. 

Algebraic diagrammatic construction (ADC) is a 
sophisticated, systematic approximation scheme for 
Green's functions. 19-21 The n-th order ADC scheme, 
ADC(n) for short, contains infinite summations of those 
classes of Feynman diagrams that derive from the first n 
orders of the Feynman-Dyson perturbation series. The 
problem of finding the pole positions and pole strengths 
of a Green's function is formulated in terms of the solu- 
tion of a Hcrmitian matrix eigenvalue problem. 

The one-particle ADC(3) scheme we use in this work 
to calculate ionization spectra employs the Dyson equa- 
tion, 13,20 which allows one to sum many diagram classes 
in addition to the summation already carried out by the 
ADC scheme itself. Exploiting the validity of the Dyson 
equation necessitates the combined treatment of the ion- 
ization potentials and the electron affinities. The latter 
are also given by pole positions of the one-particle prop- 
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agator. This approach enlarges the aforementioned ADC 
eigenvalue problem considerably. Therefore, the affinity 
block in the ADC matrix is reduced in practice by per- 
forming on this block a few (typically ten) block Lanc- 
zos iterations. 22 The most important spectral features of 
the affinity block are preserved in this way, and the fi- 
nal eigenvalue problem to be solved becomes much more 
manageable. 

The ionization potentials associated with the Xe Ad or- 
bitals in the xenon fluorides lie far above the molecu- 
lar double ionization threshold. Hence, not only two- 
hole-one-particle configurations but also three-hole-two- 
particle configurations are expected to have impact on 
the description of the core-level ionization spectra. 23-25 
The ADC(3) scheme does not contain three-hole-two- 
particle configurations explicitly, but the next-order 
scheme, ADC(4), does. Thus, one expects that the 
ADC(4) scheme yields an appreciable improvement of 
calculated core-level ionization spectra. 

The inclusion of three-hole-two-particle configurations 
enlarges the configuration space — and the associated 
eigenvalue problem — extremely. The available computer 
resources would not be sufficient to calculate ionization 
spectra of molecules like XeF 6 in an acceptable amount 
of time. The existing ADC (4) programs 23 ' 24,26 therefore 
utilize the core-valence separation approximation 25 for 
core-level ionization, which reduces the size of the con- 
figuration space. This approximation, however, makes 
it impossible to calculate ionization potentials in the va- 
lence regime. Furthermore, since core-valence separa- 
tion implies that all elements of the configuration space 
carry a core hole, those two-hole-one-particle configura- 
tions that are needed to describe core-hole decay are not 
included. For these reasons, the ADC(4) programs are 
useless for the questions we would like to address here. 
Hence, all ionization spectra were calculated within the 
ADC(3) scheme. 20 - 21 



B. Geometries and Basis Sets 

The ground-state geometries of F 2 , XeF 2 , XeF4, and 
XeFg are taken from the literature. The distance of 
the fluorine atoms in the fluorine molecule is r(F- 
F) = 1.42 A. 1 Xenon(II)-fluoride is a linear molecule 
of -Doo/i symmetry with a Xe-F distance of 1.977 A. 1 
Xenon(IV)-fluoride is square-planar (D^), the Xe-F dis- 
tance is 1.94 A. 1 All distances given refer to molecules in 
the gas phase. 

The ground-state geometry of xenon(VI)-fluoride in 
the gas phase is a distorted octahedron. 6 ' 10 It can be 
described in C% v symmetry. Nevertheless, some of our 
computations on XeF 6 were performed in Oh symmetry 
to elucidate the impact of the reduced symmetry The no- 
tation of Fig. la in Rcf. 10 is used here. Our calculations 
on XeF 6 in Cs v symmetry use the bond lengths r(Xe- 
Fl) = 1.856 A and r(Xc-F4) = 1.972 A and bond an- 
gles a — 80.8° and (i — 112.8°. These values, which were 



determined using MP2 and an all-electron basis with / 
functions, were taken from Tables 2 and 3 in Ref. 10. Our 
computations on XeF@ in Oh symmetry were performed 
at the bond length r(Xc F) = 1.952 A, which is taken 
from Table 1 in Ref. 10. 

The ADC programs 22 ' 26 (see Sec. II A) require molecu- 
lar integrals and orbital energies. These were obtained in 
Hartree-Fock calculations using the Gamess-UK 27 pro- 
gram package. The employed software does not exploit 
symmetry for the xenon atom. D 2 h symmetry is used 
for F 2 , XCF2, and XeF4. The calculation of the ioniza- 
tion potentials of XeFg in the ground-state geometry of 
C^v symmetry is performed in C s symmetry; for the ge- 
ometry of Oh symmetry, Gamess-UK assumes D2/1 sym- 
metry. Comparing the orbital energies from both the C$ v 
and the Oh geometry shows that the overall positions of 
the xenon-like orbitals are in good agreement in both 
symmetries, but the splitting of the fluorine-like orbitals 
is enlarged due to the increased interaction in Cs v sym- 
metry. 28 

The xenon and the fluorine atoms are represented by 
the DZVP (DFT orbital) 29 ' 30 basis set. The quality of 
the basis set can be estimated from Table I by comparing 
the numerically exact Hartree-Fock orbital energies (see 
the ensuing Sec. II C) of xenon with those obtained using 
the DZVP (DFT orbital) basis set. The shift of the or- 
bitals, due to the approximation introduced by the finite 
basis set, is Ae B s := £hf, numeric - £hf, dzvp- The shift 
is « 0.13 eV for the Xe4d orbitals and even less for the 
valence orbitals because basis sets are usually optimized 
with respect to the latter orbitals. This shift is neglected 
in the following, for other sources of inaccuracies are more 
significant. 



C. Relativistic Effects 

There are three main relativistic effects one has to 
take into account when examining heavy atoms like 
xenon: 31 ' 32 

1 . the relativistic radial contraction and energetic sta- 
bilization of the s and p shells, 

2. the spin-orbit splitting, 

3. the relativistic radial expansion and the energetic 
destabilization of the (outer) d and / shells. 

Effects (1) and (3) are termed scalar relativistic effects. 10 
As the theory of Sec. II A for the calculation of ion- 
ization potentials is strictly nonrelativistic, one has to 
take into account relativistic effects by a "rule of thumb" . 
This is done by performing for the xenon atom Hartree- 
Fock calculations with MCHF84 33 and Dirac-Fock calcu- 
lations, the relativistic counterpart to Hartree-Focki, 31 ' 32 
with GRASP92. 34 ' 35 Due to the spherical symmetry of 
atoms, the equations can be solved without fixed ba- 
sis sets to arbitrary precision. This gives exact sets of 
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relativistic and nonrelativistic orbitals in the mean-field 
approximation. 

In Table I, the orbital energies determined in this way 
are listed together with Hartree-Fock orbital energies ob- 
tained by a calculation employing the DZVP (DFT Or- 
bital) basis set. By comparing the orbital energies of the 
two numerical solutions of the Hartree-Fock and Dirac- 
Fock equations, one can determine the size of relativistic 
effects and correct for them in nonrelativistic computa- 
tions of the ionization potentials of the xenon fluorides. 

In order to carry out this comparison, one has to note 
that the total angular momentum j results from cou- 
pling the orbital angular momentum I with the electron 
spin: j = l±\. A Dirac-Fock calculation yields, for I > 1, 
two orbitals, 32 one for j + = 1+ | and one for j_ = I — 5. 
The j- orbital has a lower orbital energy than the j+ or- 
bital. 

The spin-orbit splitting between the j + orbital and the 
j_ orbital can be effectively removed by calculating a 
weighted mean 



£df 



(2j+ + 1) e DF ,+ + (2j_ + l)g DF ,- 
2j+ + 1 + 2j_ + 1 



(1) 



of the two Dirac-Fock orbital energies. 32 This procedure 
facilitates comparison between the Dirac-Fock and the 
Hartree-Fock calculation. The scalar relativistic shift of 
the nonrelativistic orbitals is Ae :— £df — £hf- The 
shifts are 3.077 eV for the Xe4cJ, -1.789 eV for the Xe5s, 
and -0.00283 eV for the Xe5p orbitals. The sign of the 
shift Ae is predicted by the above-given rules (1) and 
(3), because the energetic stabilization of s and p or- 
bitals leads to a lowering of the Xe 5s and Xe hp orbital 
energies and thus Ae < 0. Similarly, Ae > holds for 
the shift of the Xe Ad orbital energies due to the energetic 
destabilization of the (outer) d orbitals. 

Koopmans' theorem 15,36 states that the ionization po- 
tentials of an atom or a molecule are approximately given 
by the negative of the orbital energies, IP« — chf- As the 
relativistic shift Ae represents a relativistic correction to 
the Hartree-Fock orbital energies enF.corr = £hf+ Ae, the 
ionization potentials are corrected as follows: IP CO rr = 
IPhf - Ae. 

The mentioned shifts of the orbital energies of the 
xenon atom are used to correct the Xe4d orbital ener- 
gies for the scalar relativistic effect by adding — Ae to 
the Xe4(i ionization potential of the nonrelativistic cal- 
culations. This procedure is justified by the observation 
that the molecular Xe Ad orbitals are highly localized and 
are very similar to atomic Xe Ad orbitals. The Xe Ad ion- 
ization potentials, calculated with the ADC(3) scheme, 
take electron correlation into account. Therefore, con- 
tributions of many orbitals mix into the description of 
the Xe4d ionized states. Since the Xe4rf ionized states 
are predominantly described by the Xe4(i orbitals, the 
relativistic correction is applicable in this case as well. 

The spin-orbit splitting in the xenon atom amounts 
to 2.111 eV for the Ad orbitals and 1.436 eV for hp orbitals 
(see Table I). The spin-orbit splitting in the xenon atom 



is probably only a good approximation 6 for the Xe Ad or- 
bitals in the xenon fluorides, because the Xe hp orbitals 
suffer from a significant modification by the molecular 
bond to the fluorine atoms. The experimental value for 
the spin-orbit splitting of the Ad orbitals in the xenon 
atom is 1.984 ± 0.014eV 6 (see Table IV). The theoreti- 
cally and experimentally determined values are in satis- 
factory agreement. A good agreement is not expected be- 
cause the Dirac-Fock equations are based on a mean-field 
approximation. 32 Spin-orbit splitting is not accounted for 
in the nonrelativistic theory used throughout and is not 
considered any further in the ensuing sections. 



III. RESULTS AND DISCUSSION 
A. Mean-Field Model 

Single ionization spectra of the xenon fluorides can 
be obtained in the Hartree-Fock approximation with 
the help of Koopmans' theorem. 15,36 In this mean-field 
model, the correlation between the electrons is neglected, 
resulting in very simple spectra (see Fig. 1) that can 
provide hints for the interpretation of the more com- 
plex spectra that include electron correlation (see Fig. 2). 
The assignment of the lines in Fig. 1 to atomic orbitals 
of xenon or fluorine is not strict in the valence region, 
due to the molecular bond. In the deeper lying molecu- 
lar orbitals, this assignment is well defined. Since spin- 
orbit coupling is neglected, artificial degeneracies are in- 
troduced in the spectra of Fig. 1. The effect of the ligand 
field caused by the fluorine atoms is of course included. 

The effect of adding fluorine atoms to xenon is studied 
in Table II by Mulliken and Lowdin population analy- 
ses. 15 ' 37 Such population analyses may exhibit some ba- 
sis set dependence, 15 but the consistency between the 
two sets of results indicates reliability. One sees imme- 
diately that charge is withdrawn from the xenon atom 
to the fluorine atoms: 1.1 electron charges in XeF 2 , 1.9 
in XCF4, and 2.3 in XeF 6 . Due to the Cs v symmetry of 
the ground-state geometry of XeF 6 , there are two kinds 
of fluorine atoms, with different distances to the central 
xenon atom. These differences in nuclear separation are 
reflected in Table II by the fact that fluorine atoms which 
are further away from the xenon atom acquire less charge. 

At first sight one might assume that the computed 
change in charge density involves only the valence elec- 
trons and has little effect on the inner molecular orbitals 
of the xenon fluorides. In the hydrogen atom, the wave 
functions of the higher lying shells, which are unoccupied 
in the ground state, have a considerable amplitude in the 
spatial regions of the lower lying shells of the same angu- 
lar momentum. 38 As this argument also holds in the case 
of the xenon fluorides, a reduction of valence electron 
density on the xenon atom leads to a less efficient screen- 
ing of its nuclear charge and consequently to a lowering 
in energy of the energetically low lying molecular orbitals 
with a dominant contribution on the xenon atom. 
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tttv _ ; ^ „ 1 

nr orbital 


£hf, dzvp 


^HF, numeric 


Dr orbital 


£DF 


£df 


4d 


-2.78280 


-2.77788 


4d 3 / 2 
4d 5 / 2 


-2.71133 
-2.63376 


-2.66479 


5s 


-0.946253 


-0.944414 


5s 


-1.01014 


-1.01014 


5p 


-0.457894 


-0.457290 


5P1/2 
5P3/2 


-0.492572 
-0.439805 


-0.457394 



TABLE I: Hartree-Fock (HF) and Dirac-Fock (DF) orbital energies of xenon. Hartree-Fock orbital energies are given for the 
P (DFT 

' : 35 £df is calculated using Eq. (1). All data are given in Hartree. 



DZVP (DFT orbital) 29 ' 30 basis set and for the numerically exact solution. 33 The Dirac-Fock orbital energies also are numerically 



exact. 
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FIG. 1: (Color) Single ionization spectra of Xe, F2, XeF2, XeF4, and XeF6 calculated at the Hartree-Fock level (Koopmans' 
theorem) . The assignment of the lines to atomic orbitals of xenon or fluorine origin in XeF„ is not strict in the valence region, 
due to the molecular bond. 



Conversely, the increase of valence electron density on 
the fluorine atoms leads to a more efficient screening of 
their nuclear charge and consequently raises the energy of 
the lower lying molecular orbitals with a dominant con- 
tribution on fluorine atoms. The charge withdrawn from 
the xenon atom is shared among several fluorine atoms. 



The net increase of charge density on each fluorine atom 
is w 0.5 electron charges. This value is, of course, much 
smaller than the loss of charge density on the xenon atom 
and hence the impact on the inner molecular orbitals of 
fluorine character is clearly much smaller than for xenon. 
It is most dramatic in XeF2 where each fluorine atom ac- 
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Ae 
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12.10 


XeF 2 


12.35 


12.76 


XeF 4 


13.1 


13.07 


XeF 6 


12.35 


12.56 



TABLE II: Mulliken and Lowdin population analyses of XeF2, 
XeF 4 , and XeF 6 . 



quires the largest fraction of charge. The screening of the 
nuclear charge is largest in XeF 2 and decreases in XeF 4 
and XeF 6 . 

The effects of this model can be seen in Fig. 1. The 
positions of the Xe 5s and Xe Ad lines shift to higher bind- 
ing energies with an increasing number of fluorine atoms. 
The energy differences between the corresponding lines 
of XeF 2 and XeF 4 and of XcF 4 and XeFg are nearly 
equally large. The F 2p and F 2s lines also shift slightly 
to higher binding energies with increasing number of flu- 
orine atoms. The charge redistribution over several flu- 
orine ligands we mentioned earlier furnishes an explana- 
tion for the shift of the fluorine lines. This interpretation 
is further supported by a comparison with F 2 . The mean 
of the F 2p lines in F 2 and the mean of the F 2s lines in F 2 
arc higher in energy than the mean values of the corre- 
sponding lines in XcF 2 . The first ionization potential is 
nearly constant for all xenon fluorides studied. Its value 
is « 12.5 eV. 

In the fluorine molecule, the two F2s lines are split 
considerably due to the molecular bond. The splitting 
of the two F 2s lines in XcF 2 is tiny due to the large 
separation of the fluorine atoms. In XeF 4 and XeF 6 , the 
fluorine atoms are closer and interact. This results in 
a larger splitting of the F2s lines in comparison to the 
splitting in XeF 2 . 

The Xe Ad lines are quintuply degenerate in the single 
ionization spectrum of the xenon atom in Fig. 1. In XeF 2 , 
the degeneracy is lifted by the ligand field of the fluo- 
rine atoms and three distinct lines become visible. The 
three lines reflect the spatial orientation of the Ad or- 
bitals. There is a ligand field along the molecular axis. 
Perpendicular to the molecular axis there is no shift, re- 
sulting in a total of three lines. In XeF 4 there is only 
one dimension left that is unaffected by the ligand field: 
the axis perpendicular to the molecular plane. The spec- 
trum, Fig. 1, shows that there are two distinct groups of 
Xe4cZ lines in this case. 

In XeFg, the situation changes because the fluorine lig- 
ands are grouped around the xenon atom in such a way 
that XeF 6 is close to octahedral symmetry. Oh symme- 
try in XeFg leads to two distinct Xe4<i orbital energies 
separated by only 0.1568 eV. Lowering the symmetry of 



TABLE III: Comparison of the calculated lowest (first) ion- 
ization potentials (IP) of Xe, XeF2, XeF 4 , and XeF6 with 
experimental results. The first ionization potential of xenon 
is taken from Ref. 1. The other first ionization potentials are 
taken from Ref. 4. All data are given in electronvolt. 



XeF 6 to Czv results in three distinct XeAd orbital ener- 
gies with a total splitting of 0.3682 eV, which is still small 
compared to the much larger Xe Ad splittings in XeF 2 and 
XeF 4 . 



B. Correlation Effects 

Going beyond the Hartree-Fock description of the 
molecules by using ADC(3) leads to the ionization poten- 
tials plotted in Fig. 2. Figure 1 helps identify the one- 
particle origin of the states in Fig. 2. The Xe4c? lines 
are energetically clearly separated from the outer and 
inner valence lines in all spectra. They are located be- 
tween 72 eV and 82 eV. The F 2s-derived lines are present 
between 35 eV and 45 eV. Correlation effects do not 
change the fact that the lowest ionization potentials are 
approximately equal in all compounds. 

In fact, at the lowest ionization potentials, which are 
associated with the Xe5p levels, there is qualitative 
agreement between the correlated and the mean-field cal- 
culation: The pole strengths are very close to unity; 
the one-particle picture is valid. Nevertheless, ADC(3) 
improves on the Hartree-Fock method in a quantitative 
manner. In the correlated calculation, the energetically 
lowest cations are stabilized by about 1 eV relative to the 
respective one-particle energies. The first ionization po- 
tentials of xenon and its fluorides, as calculated within 
the ADC(3) scheme, are compared to the experimental 
results in Table III. The agreement of experimental and 
calculated ionization potentials is good. 

At higher energies, there is typically not even qualita- 
tive agreement between mean-field and many-body treat- 
ment. Overall, there is a dramatic change of the ioniza- 
tion spectra due to the inclusion of many-body effects. 
The importance of electron correlation grows with in- 
creasing number of fluorine atoms. This is reflected in 
energy shifts and reduction of intensity in the outer va- 
lence part of the spectra; in the appearance of numerous 
satellite lines in the vicinity of 20 eV and above; and, in 
particular, in the eye-catching breakdown of the molec- 
ular orbital picture of ionization seen in the Xe 5s and 
F 2s inner-valence region. 

The latter phenomenon is common in molecules: 39 The 
intensity originally confined to a single orbital in Fig. 1 
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FIG. 2: (Color) Single ionization spectra of Xe, F2, XeF2, XeF4, and XeF6 calculated with the one-particle ADC(3) program 
(see the text). The assignment of lines to atomic orbitals of xenon or fluorine origin is not done for the XeF„ molecules. 



is spread over many cationic states (see Fig. 2). The 
breakdown of the molecular orbital picture of ioniza- 
tion is caused by compact two-hole-one-particle config- 
urations that are close in energy to a one-hole state. 39 
The strong Coulomb coupling of these configurations to 
one-hole states leads to a broad shape of lines, which 
has been discovered in the inner-valence region of many 
molecules. 39 

Due to the many-body description of a molecule by 
the ADC scheme, decay processes of electronic states of 
the ion that lie above the double ionization threshold 
also become describable. 40 In particular, the final state 
of the electronic decay of an ionized molecule can be ap- 
proximated in terms of two-hole-one-particle configura- 
tions describing the emitted electron (the particle) and 
the residual dication (the two holes). 41 The shape of each 
decaying state can be identified as a thin bundle of lines, 
which mimic a discretized Lorentzian curve. 40 



In order to unambiguously identify a decaying state, 
the breakdown of the molecular orbital picture of ioniza- 
tion has to be separated clearly from the decaying states 
in a small energy range where each decaying state is rep- 
resented by a discretized Lorentzian curve. Furthermore, 
there can be a mixing of both phenomena for a certain 
state, but it is hard to identify a decay curve for one-hole 
states that suffer from the breakdown phenomenon. 

The inner-valence region of the xenon fluorides com- 
prises the Xe 5s and F 2s states. The F 2s states are 
subject to the breakdown phenomenon and can also de- 
cay by emitting an electron. 28,42 The Xe 5s one-hole state 
suffers solely from breakdown. No decay is possible be- 
cause its ionization potential is below the double ioniza- 
tion threshold. 28,42 Many two-hole-one-particle configu- 
rations can couple to the Xe 5s one-hole state (see Fig. 2). 
As the Xe 5s line is close to the double ionization thresh- 
old, 28,42 many two-hole-one-particle configurations are 



7 



close in energy to this state: The excited states are gen- 
erally very dense in this energetic region. 

The spectra exhibit a considerable increase of break- 
down of the molecular orbital picture for the inner va- 
lence with an increasing number of fluorine atoms. Ob- 
viously, the number of suitable configurations rises due to 
the addition of fluorine atoms. Since xenon fluorides are 
very symmetric, these fluorine atoms are all equivalent 
(except for XeFg in Cs v symmetry), and the number of 
equivalent two-hole-one-particle configurations with one 
hole on the xenon atom and one hole on a fluorine atom 
doubles between XeF 2 and XeF 4 . In XcF 6 there are even 
more such configurations than there are in XcF 4 , but 
they are no longer equivalent due to the Cs v symmetry 
of the XeFg ground state. 

Another contribution to the increase of one-hole states 
involved in the breakdown phenomenon can be attributed 
to the decrease in symmetry: spherical symmetry (Xe), 
Dooh (XcF 2 ), D ih (XeF 4 ), C 3v (XeF 6 ). This decrease 
in symmetry leads to an increase in excited electronic 
configurations that can couple to the respective one-hole 
states. 

The one-particle picture provides a rather good de- 
scription for the Xe4d lines. Two-hole-one-particle con- 
figurations close to the energy of the core one-hole config- 
urations are spatially diffuse and thus couple only weakly 
to the core one-hole configurations. Therefore, in the 
core regime, the mixing of two-hole-one-particle config- 
urations with one-hole configurations is relatively weak. 
The splittings of the Xe Ad lines in the xenon fluorides are 
of comparable size to those in Fig. 1, and the individual 
lines can be identified easily in all molecules but XeFg. 

Since in all xenon fluorides the Xe4<i lines are above 
the double ionization thresholdi, 28,42 Xe 4e?-ionized xenon 
fluorides can undergo Auger decay. The Xe4c? lines in 
Fig. 2 show signs of decay by exhibiting features of dis- 
cretized Lorentzian curves. The shape of the approxi- 
mate Lorentzian curve representing a decaying state in 
the spectra depends highly on the number of two-hole- 
one-particle configurations in the energy range of the de- 
cay electron. 

The decay electrons that result from the decay of a 
Xe Ad vacancy are highly energetic (up to 60 eV) . 
Hence, a satisfactory description in terms of Gaussian 
basis sets is not feasible. For this reason, the Xe4d ion- 
ization potentials in Fig. 2 correspond in our description 
essentially to one-hole states (a pole strength close to 
unity). Nevertheless, the number of states which orig- 
inate from two-hole-one-particle configurations around 
the XeAd ionization potentials grows substantially with 
the number of fluorine atoms. In XeFg one observes very 
densely lying states around 83 eV with a high contri- 
bution of two-hole-one-particle configurations. The ex- 
treme change in the importance of two-hole-one-particlc 
configurations in XeF 6 compared to XeF 4 , XeF 2 was in- 
vestigated by us. 

The six fluorine atoms in XeFg form a distorted oc- 
tahedron around the xenon atom. The basis set of the 



fluorine atoms may therefore help improve the descrip- 
tion of the decay electron from Xe Ad- ionized XeFg con- 
siderably. To test this assumption, we calculated the ion- 
ization spectrum of XeF 2 with the DZVP (DFT orbital) 
basis set augmented by a few diffuse functions. The re- 
sulting spectrum did not change much compared to that 
of Fig. 2. 

To test the effect of the basis functions on the fluorine 
atoms, we performed another calculation of the ionization 
spectrum of XeF 2 with fluorine basis functions attached 
to ghost centers arranged such as to yield an octahedron 
of fluorine basis functions surrounding the central xenon 
atom. The extra basis functions had a minor effect on 
the spectrum and did not account for the drastic effect 
observed in Fig. 2 for XeFg. 

These investigations convinced us of the suitability 
of the DZVP (DFT orbital) basis set for our concerns 
but did not clarify the reasons for the drastic change in 
Fig. 2. To this end, we calculated the ionization spec- 
trum of XeFg in Oh geometry. The resulting Xe 4c? lines 
had a shape similar to that of the Xe4c? lines in XcF 4 , 
revealing the great importance of the decreased symme- 
try for the description of XeFg in our computations. The 
reduction in symmetry from Oh to Cs v enables a vastly 
enlarged amount of two-hole-one-particle configurations 
to couple to the Xe4<i one- hole configurations. To some 
extent, this is similar to an increased basis set, as a mul- 
titude of configurations are provided additionally for the 
description of the ionized states. For sure, this improves 
the description of those decay electrons with low kinetic 
energy. 

The charge transfer from the xenon atom to the fluo- 
rine atoms in the ground state, discussed in Sec. Ill A, 
leaves a positively charged central xenon atom. There- 
fore, it is not surprising that XeF 4 and XeFg possess 
positive electron affinities, i.e., (XeF 4 )~ and (XeFg) - 
are electronically stable anions. We have computed the 
respective electron affinities by ADC(3). The value ob- 
tained for XeF 4 (a doubly degenerate state) is 0.66 eV; 
the electron affinities of XeFg lie at 1.30 eV, 1.31 eV, and 
1.81 eV, respectively. The figures are presumably not 
very accurate because the DZVP (DFT orbital) basis set 
is not particularly well suited for describing the diffuse 
states of anions. In XeFg, the fluorine atoms form a 
cage of negative charge, of nearly octahedral symmetry, 
surrounding the central, positively charged xenon atom. 
The resulting potential well is similar to the model po- 
tential well discussed in Refs. 28,43. In these references, 
both bound states as well as scattering resonances in the 
model potential well were investigated. 

The Xc4g? and F2s lines in Fig. 2 are shifted to 
higher ionization potentials with increasing number of 
fluorine atoms (see Sec. Ill A). This general trend, seen 
in the mean-field approximation of Fig. 1, can still be 
identified if electron correlation is taken into account. 
Nevertheless, the difference between the shifts of the 
Xe4cZ lines of XeF 2 and XeF 4 is larger than that be- 
tween the Xe Ad lines of XeF 4 and XeFg . Electron cor- 
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Cmpd. 


Line 


IPcxpt 


IPadc 


IPrel 


Tcxpt 


Xe 


1 


69.525 (10) 


72.90 


69.82 


0.207 (4) 




2 


67.541 (9) 






0.202 (4) 


XeF 2 


1 


72.568 (6) 


76.52 


73.44 


0.248 (8) 




2 


72.248 (6) 


76.14 


73.06 


0.223 (10) 




3 


70.601 (13) 


75.63 


72.55 


0.264 (26) 




4 


70.421 (9) 






0.256 (27) 




5 


70.179 (6) 






0.214 (19) 


XeF 4 


1 


75.098 (6) 


79.76 


76.68 


0.319 (8) 




2 


74.729 (7) 


79.60 


76.52 


0.255 (8) 




3 


73.140 (10) 


79.01 


75.93 


0.392 (10) 




4 


72.816 (10) 


78.89 


75.81 


0.210 (27) 




5 


72.661 (5) 






0.225 (26) 


XeF 6 


1 


77.462 (13) 


80.86 


77.78 


0.32 (4) 




2 


77.321 (11) 


80.59 


77.51 


0.25 (3) 




3 


75.53 






0.33 




4 


75.38 






0.25 




5 


75.25 






0.25 



TABLE IV: Peak positions and widths of the Xe4d lines 
in Xe, XeF2, XeF.4, and XeF6. Labels — Cmpd: compound; 
IPexpt: experimental peak position; IPadc: calculated peak 
position; IP rc i: calculated peak position with relativistic cor- 
rection; r expt : experimental peak width. IP CX pt and r oxp t are 
photoelectron-experimental data, reproduced from Table 1 in 
Ref. 6. For XeF6, the data with an experimental resolu- 
tion of 0.11 eV are taken. The values in parentheses are the 
standard deviations of the peak positions and peak widths. 6 
IPadc and IP rc i are sorted descending in energy for each com- 
pound. The value of IP exp t does not necessarily correspond 
to the values of IPadc and IP ro i in the same row as spin-orbit 
splitting is neglected to obtain the latter ones. All data are 
given in electronvolt. 



relation reduces the effect caused by charge depletion on 
the xenon atom. 

The nonrelativistic IPadc are listed in Table IV to- 
gether with the IPrel obtained by applying to IPadc the 
relativistic correction — Ae of Sec. II C. The number of 
distinct Xe Ad lines is smaller in the nonrelativistic spec- 
tra due to a higher degeneracy caused by neglecting the 
spin-orbit coupling. The Xe Ad lines of XeFg cannot be 
identified clearly in Fig. 2. One expects from Sec. Ill A 
two lines for a ground-state geometry of Oh symmetry 
and three lines for a ground-state geometry of C% v sym- 
metry. The two ionization potentials, in the energy range 
of the XeAd lines of XeF6 (80-83 eV), with maximum 
pole strength are listed in Table IV. 

The Xe Ad IP ro i differ from the experimentally obtained 
data in Table IV by only 1.5-2.5 eV (w 3%), which is a 
good agreement in view of the complexity of the problem 
and the necessity to describe the decay simultaneously. 
The reason for the deviation is twofold. First, the spin- 
orbit splitting is neglected, which amounts to 2.111 eV 
for the XeAd lines of the xenon atom, according to 



Sec. II C. Second, since the calculations are based on the 
ADC(3) approximation, the three-hole-two-particle con- 
figurations are neglected [ADC (4) calculations are cur- 
rently beyond reach]. The inclusion of these additional 
configurations would shift the ionization potentials of the 
XeAd lines further to lower energy, due to an improved 
description of core-hole relaxation. 23 

Let us now turn our attention to the last column in 
Table IV: the experimental widths of the XeAd photo- 
electron peaks. 6 It is evident that, on average, the widths 
increase from Xe to XeF 6 . According to Cutler et al., 6 
this effect probably cannot be explained in terms of vi- 
brational broadening. If the data reflect more or less 
pure Auger widths, then the experimental findings are 
in contradiction to expectations based on the one-center 
model. 44 The electronegative fluorine ligands withdraw 
charge from xenon (see Table II) . Correspondingly, fewer 
valence electrons are available for a local, atomic mech- 
anism of Xe Ad Auger decay. Hartmann's semiempirical 
multi-center model, 45 on the other hand, predicts a larger 
C Is Auger width for CF 4 than for CH4, which, though 
seemingly in line with the trend in XeF„, is inconsistent 
with experiment. 46 

In our work on neon clusters, we demonstrated not 
only that a Ne 2 s hole in a cluster can decay by electron 
emission 47 — a fact confirmed by a recent experiment 48 — 
but also that the decay lifetime depends sensitively on the 
number of nearest neighbors surrounding the atom carry- 
ing the inner- valence hole. 41 The Ne 2s lifetime drops in 
a monotonic manner from about 80 fs in NeJ to less than 
5 fs in Nejj. The mechanism responsible for this dramatic 
effect is referred to as interatomic Coulombic decay. In 
heteroclusters, additional electronic mechanisms play a 
role. 49 We will demonstrate in an upcoming paper 42 how 
the decay mechanisms discovered in clusters elucidate the 
counterintuitive behavior of the Xc4d Auger widths in 
XeF„. 

Here, we would like to emphasize that traces of a purely 
electronic origin can already be found in our ADC (3) 
spectra, Fig. 2. As the number of fluorine ligands grows, 
there are more and more two-hole-one-particle configura- 
tions that can couple to the Xe Ad holes. In other words, 
the number of dicationic decay channels grows, and there- 
fore we observe in Fig. 2 an increasing loss of main-line 
intensity and the emergence of dense line bundles. This 
is the signature in our finite-basis-set calculations of en- 
hanced electronic decay, 40 as pointed out previously. 



IV. CONCLUSION 

The xenon fluorides are very special systems. Not only 
are they rather peculiar chemical compounds, they are 
also particularly well suited to investigating the system- 
atic evolution of many-body effects. 

While outer-valence ionization is mostly compatible 
with the mean-field approach, the situation is strikingly 
different in the inner-valence region. At the Hartree-Fock 
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level, the chemical environment leads to mere shifts of the 
Xe5s and F2s lines. The ADC(3) calculation, however, 
reveals that electron correlation can induce additional, 
pronounced environmental effects. Xe 5s ionization in 
XeF2 can be relatively well understood in terms of a one- 
particle picture. In XCF4, this picture begins to crumble, 
and in XeF 6 , there is not a single cationic eigenstate that 
overlaps well with the Xc 5s one-hole configuration. 

The F 2s states are even more interesting, since they 
are higher in energy than the double ionization thresh- 
old. In the xenon fluorides, electron correlation enforces 
the excitation of a large, system-dependent number of 
F 2s-derived electronic states. All of these states are 
resonances decaying by electron emission. Neither the 
decay mechanism nor the decay time scale are known. 
Resolving these issues may have important consequences 
for other molecules. This poses an exciting challenge to 
both theory and experiment. 

This is the first paper dealing with the inner-valence 
physics of XcF„ at a correlated level. We are not aware 
of any experimental spectra in this energy region. The 
quality of our calculations can be assessed, however, by 
reference to our ADC (3) results for the Xe4d core lines. 
Upon correction for scalar relativistic effects, our numeri- 
cal data agree with experiment within a few percent. The 



remaining difference can be attributed to an incomplete 
description of core-hole relaxation and to our neglect of 
spin-orbit coupling. The ADC(3) scheme is capable of 
reproducing the experimentally observed line shifts as 
well as the counterintuitive Auger broadening. A detailed 
analysis of the latter phenomenon will be presented else- 
where. 
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